library(tidyverse)
library(sf)
library(tigris)
library(censusapi)
library(leaflet)
Sys.setenv(CENSUS_KEY="b88495f8315f45b2d100dee9ba8f4c489a7371c2")
blockgroups <- block_groups("CA",037,progress_bar = F)
# list <- c("060372060102","060371997001","060371997002","060379800101","060371977003","060372071032","060372060101")
#
# narrow_blockgroups <-
# blockgroups %>%
# filter(GEOID %in% list)
test <- getCensus(
name = "acs/acs5",
vintage = 2018,
region = "block group:*",
regionin = paste0("state:06+county:037"),
vars = "B01003_001E"
) %>%
transmute(
origin_census_block_group = paste0(
state,county,tract,block_group
),
pop = B01003_001E
)
test1 <- test %>% left_join(blockgroups,by = c("origin_census_block_group" = "GEOID")) %>% dplyr::select(origin_census_block_group, pop,geometry) %>% st_as_sf()
leaflet() %>%
addTiles() %>%
addPolygons(
data = test1,
fillColor = "blue",
color = "white",
opacity = 0.5,
fillOpacity = 0.5,
weight = 1,
label = ~pop,
highlightOptions = highlightOptions(
weight = 2,
opacity = 1
)
)